Mechanism of Ferroelectricity in d 3 Perovskites - a Model Study 
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By means of a model Hamiltonian approach we study the role of volume expansion, Hund's cou- 
pling and electron correlation in the standard hybridization mechanism for ferroelectricity in cubic 
CaMn03, a prototypical non-d° perovskite. Our results establish that the ferroelectric instability 
arises from a subtle balance between different energy contributions, explaining the origin of its en- 
hancement under negative pressure. Expansion of volume is found to cause a strong reduction of 
the elastic energy, while leaving almost unchanged the tendency of Mn states to form covalent bonds 
with the surrounding oxygens. Hund's coupling with local spins of magnetic cations can reduce and 
even suppress the instability towards the ferroelectric state. 
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I. INTRODUCTION 

Since the revival of interest in the magnctoclcctric ef- 
fect due to the possibility of multiferrocitj** with simulta- 
neous presence of ferroelectric and magnetic order, many 
attempts have been made to create multiferroic materi- 
als focusing on the wide variety of transition metal ABO3 
perovskites. 

The classical ferroelectrics as well as many magnetic 
materials belong 2 ^ to this class of compounds. However, 
overlap between ferroelectric perovskites and magnetic 
perovskites are rare, with the exception of compounds 
like BiFe0 3 and B1M11O3. This lead to the suggestion 
of general "exclusion" rule for perovskites^—, according 
to which ferroelectricity occurs for d° B site transition- 
metal (TM) cation^, while the magnetism requires d n B 
site cations with n ^ 0. This mutual exclusion of fer- 
roelectric and magnetic properties is rationalized by the 
cases of conventional ferroelectric perovskites like, e.g., 
BaTi0 3 , KNb0 3 or Pb(ZrTi)0 3 , for which the driving 
force behind the ferroelectric distortion comes from the 
tendency of TM empty d states to establish a strong co- 
valency with the surrounding oxygensi^ The hybridiza- 
tion driven mechanism is maximally efficient when the 
antibonding states, with predominant TM d character, 
arc completely empty, hence for d° B site perovskites. 
The above mentioned exclusion principle is violated for 
cases like BiFe0 3 and BiMn0 3 , for which the polarization 
arises in A sublatticc due to the lone-pair stereochemical 
activity of Bi 6s electrons^ and magnetism arises in B 
sublattice, with no apparent coupling between the two or- 
der parameters. For technological purpose, on the other 
hand, it is highly desirable to invoke coupling between 
the two order parameters, so that magnetism can be con- 
trolled by electric field and polarization by magnetic field. 
Along this line of thought, a number of magnetically in- 
duced multiferroic perovskites have been synthesized and 
studied, in which the complicated magnetic structure ei- 
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ther with non-collinear spiral spin configuratio: 
collincar ttii E-typc antifcrromagnetio^— structure, 
breaks the inversion symmetry and induces polarization. 



Also non-trivial charge distributions in systems display- 
ing both charge and/or orbital order may significantly 
contribute to the onset of polar ferroelectric phases^—, 
even concomitantly with magnetism. These compounds, 
known as improper ferroics, normally do not exhibit off 
centric movement of magnetic cations. While such elec- 
tronically driven mechanism has generated lot of aca- 
demic interest, the resulting polarization is often tiny, 
precluding technological applications. Search is, there- 
fore, on to obtain multiferroic materials showing magne- 
toelectric coupling with reasonably high transition tem- 
perature and polarization values. 

Since the conventional lattice driven mechanism of po- 
larization with off centric movement of B site cations 
arises due to delicate balance between the hybridization 
energy gain and the repulsive electrostatic energy, the 
general applicability of the "exclusion principle" for a 
d n perovskite with n < 5 is difficult to predict a-priori. 
d n perovskite, if can be made to a conventional ferro- 
electric, is expected to exhibit magnctoclcctric coupling 
as well as presumably large values of polarization. Mn 
based perovskites like CaMn0 3 , BaMn0 3 or SrMn0 3 
may be considered as possible candidates, since B site 
cation in this case is in a d 3 electronic configuration, with 
electrons occupying, within a cubic crystal field, the 3- 
fold degenerate t2 g states, leaving the e g states unoccu- 
pied. Noticeably, like conventional d° ferroelectric per- 
ovskites, they exhibit anomalously large Born effective 
charges (BECs) of Mn ionsi^Sr— The study by Ederer et 
al£2- showed that the largest anomalous contributions to 
BECs can be understood in terms of a charge transfer be- 
tween TM d and O p states along the a— type channel, 
suggesting that the hybridization mechanism is mostly 
e g -driven. Attempts have been made to induce ferro- 
electricity in these compounds through increase of vol- 
ume by applying negative or chemical pressur o 21 ' 22 , or 
in epitaxially strained thin films'. Interestingly, the po- 
larization in these predicted systems may be larger than 
that typically expected in magnetically-induced multi- 
ferroics such as TbMn0 3 (P < 0.1 ^C/cm 2 ) or HoMn0 3 
(P ~ 0.5 — 2 /iC/cm 2 ). Rather large polarization value 
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of P = 54 pC/cm 2 has been predicted for cpitaxially 
strained SrMn0 3 ^ 

A microscopic understanding of this interesting behav- 
ior demands further study, in particular what makes fer- 
roelectric instability stronger at larger lattice constants. 
Hybridization effect is expected to be weaker for in- 
creased TM - O bondlengths and therefore would dis- 
favor the hybridization-driven ferroelectricity even fur- 
ther. The role of Hund's coupling between t^g and e g 
spins, which has been suggested^ to prevent ferroelec- 
tricity in compounds like CaMn03, needs to be also un- 
derstood in the context of increased lattice constant. In 
the present study, we investigated this issue in a model 
Hamiltonian based approach. The model Hamiltonian 
approach, as opposed to the first-principles density func- 
tional theory (DFT) calculations, has the advantage of 
analysing the different energy contributions individually 
and providing a microscopic understanding of the influ- 
ence of various different factors. However, since the fer- 
roelectric instability arises due to delicate balance be- 
tween different competing energies, we employ DFT tools 
for accurate estimates of different energies. Specifically, 
we have used muffin-tin orbital (MTO) based NMTO- 
downfolding technique^ to estimate the hopping inter- 
actions between nearest neighbor Mn and oxygen atoms, 
and that between next nearest neighbor oxygens. The en- 
ergy associated with the stiffness of the off-centric, ferro- 
electric movement has been estimated using total-energy 
DFT calculations carried out employing the plane-wave 
pseudopotential method as implemented within VASPj2£ 
We assumed the high temperature cubic crystal structure 
of CaMn03, which is the prototypical cubic TM per- 
ovskitc crystal structure. We have not considered in our 
study the antiferrodistortive tilt and rotation of MnC>6 
octahedra which have been foun d 21 ' 22 to compete with 
the ferroelectric instability. Our study, therefore, focuses 
on the competition between the pair-breaking Hund's 
rule coupling, the hybridization effect and the short- 
range, electrostatic repulsion in destabilizing or stabi- 
lizing ferroelectricity in a prototypical cubic perovskitc 
structure, and the effect of volume expansion on it. 



II. METHODS AND COMPUTATIONAL 
DETAILS 

For our study, we used a minimal model designed to 
reproduce the tendency of empty e g states to establish 
covalcnt bonds with surrounding oxygens. Neglecting 
the tig electrons motion, the model should describe the 
hopping between filled oxygen p state and empty Mn 
two-fold e g states on the background of localized tig 
spins. Within the cubic crystal field, the e g states es- 
tablish prevalently a bonds with O p states parallel to 
the bond direction^, as confirmed by Wannier function 
analysi o 20 i 27 . We therefore considered a 5-band d — p 



model, given by H = Hd + H p + Hd p , where 

i.a,a i 

takes into account the local physics on TM site, 

s • X'Mwv. (2) 

describes the oxygen sites, and 

H d P = E V&dLPj* ( 3 ) 

describes the hybridization between Mn and O sites. 
d\ creates an electron with spin a on a Mn at site 
% and in the orbital state x 2 — y 2 (3z 2 — r 2 ) for a = 1(2), 
while Pj a creates an electron on O site j = «±e/2, with e 
the direction unit vector x, y, z, in the orbital state x, y, z 
respectively. £q is the d(p) onsite energy of Mn e g (O- 
p), whereas t^, and are hopping parameters that 
describe hybridization between neighboring O p states 
and between Mn e g states and O p states, respectively. 

In order to get a faithful estimate of the onsite ener- 
gies and the hopping interactions, we have carried out 
NMTO-downfolding calculations 2 ^. Starting from a full 
DFT calculation for cubic CaMn0 3 , NMTO-downfolding 
arrives at a 5-orbital Hamiltonian by integrating out de- 
grees of freedom involving Mn-^g, O-p^ and those in- 
volving Ca. The downfolded 5-orbital Hamiltonian is 
built by defining energy-selected, effective Mn-e s and O- 
p a orbitals which serve as Wannier-likc orbitals span- 
ning the downfolded five bands. Linear muffin-tin cal- 
culations (LMTO)2£ were performed in order to evaluate 
self-consistcntly the potential parameters needed for the 
NMTO procedure. 

The second and third term in Eqn. ([]} describe 
the Coulomb interactions in the 3d shell of the TM 
cation. Treating the ti g electrons as classical spins Si, 
the interaction between them and the e g electrons can 
be accounted for by the Hund's coupling Jh, connect- 
ing the localized core spins Si with the spins Si = 

Y, a ,a,a' Aaa W d i a a' of tnc mobile electrons, f being 
the Pauli matrices, as given in the second term in Eqn. 
([1]). In most of our calculations, we have assumed the an- 
tiferromagnetic G-type arrangement of ti g spins, though 
to study the influence of spin arrangements, we have also 
considered the ferromagnetic alignment of ti g spins. The 
third term, H e - e describes the Coulomb interactions be- 
tween e g electrons, as given below. We drop the site 
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index for sake of simplicity. 

H e - e = - ^ n aa n aS + (4) 

a,£T 

JJI 

C y ^ ^ (l — ^aa') footer f^a. 1 a' 
a, a' ,cr,(7' 

(1 - <W)(1 - <W) 4<T rf L' d a'<T' d a'<Tj > 

U,U',J are the intraband Coulomb, interband Coulomb 
and interband exchange couplings respectively^, defined 
within the e s manifold. The values of U, U', J were esti- 
mated as, U = 7.02 eV, U' = 4.97 eV and J = 1.025 eV, 
as relevant for CaMnOs*^ Treating H e - e at the mean- 
field level, applying a standard Hartree-Fock linearization 
of two-body operators, the problem reduces to a single- 
particle problem, which can be solved straightforwardly 
in reciprocal space. Considering the fact that the e g man- 
ifold is empty, such a mean-field treatment is justified. In 
order to take into account the G-type antiferromagnetic 
configuration of ti g spins, we defined a unit cell hosting 
two TM cations and six O ions^i. After transformation 
in reciprocal space, the Hamiltonian can be diagonalized 
for any given value of momentum fe; integration was then 
performed using a nested Clenshaw-Curtis quadrature 
rule^. 

In order to investigate the ferroelectric instability, we 
considered an off-centering displacement u of each Mn 
cation along the +z direction, i.e., towards one of the 
apical oxygen, as considered in Ref. [20l The correspond- 
ing hopping parameters were modified proportionally to 
~ (d/do)~ 3 - 5 , d and do being the distorted and undis- 
torted Mn-0 bond lengths^ respectively— The elastic 
energy loss associated with the distortion, u, was evalu- 
ated using total-energy DFT calculations, as mentioned 
earlier. The plane wave calculations, carried out for 
this purpose, were performed using projector augmented 
wave (PAW) potentials^ and the wavefunctions were ex- 
panded in the plane wave basis with a kinetic energy 
cut-off of 450 eV. Reciprocal space integration was car- 
ried out with a k-space mesh of 6x6x6. 



III. RESULTS 

A. Volume effects on covalency energy, role of 
correlation and Hund's coupling 

Using the NMTO-downfolding technique, as explained 
in Section II, we evaluated the onsite energies and various 
hopping interactions, relevant for CaMn03 for two cu- 
bic crystal structures, one at equilibrium volume^ with 
lattice constant ao = 3.73 A and another at increased 
volume, as under a negative isotropic pressure, with 2% 



increased lattice constant ao = 3.805 A. The obtained pa- 
rameters are listed in Table HI The hopping interactions 
between Mn and O sites is given by, V°- = (— l) Mi i t^, 
where .\/, , = 2 if j = i + e/2 and .\/, , = 1 if j = i - e/2; 
P = x, y, z labels the p orbital state that hybridize with 
Mn e g states. Within a cubic symmetry, the interaction 
between TM x 2 — y 2 and O p z is zero, and therefore not 
listed in the table. 
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TABLE I: Onsite and hopping parameters at equilibrium 
and 2% increased lattice constant, obtained from NMTO- 
downfolding calculations. 

Naively, one would expect that the tendency to form 
strong covalent bonds would decrease in the expanded 
volume structure. As is seen from the parameters listed 
in Table HI the hopping parameters are reduced by ~ 10% 
upon 2% increase in lattice constant. However, inter- 
estingly, the volume expansion also influences the level 
splitting between d and p states, as is seen in Table HI 
Following Khomskii^, one can estimate the energy gain 
due to the off-centering of the TM ion, as described be- 
low. When a Mn atom shifts towards an apical oxygen, 
reducing the corresponding Mn-0 bond length, it forms 
a strong covalent bond with this particular oxygen at the 
expense of weakening the bonds with other oxygen ions. 
Within a linear approximation and neglecting the bonds 
between the Mn ion and the in-plane oxygens, whose 
bond lengths are only slightly increased in the consid- 
ered off-center shift, the covalency energy gain can be 
estimated as SE cov —2 (u< p( j) 2 /A, A being the charge- 
transfer energy, given by \sd — e p \. An estimate of the 
changes in covalency energy gain induced by volume ex- 
pansion may be obtained by considering the ratio t 2 d / A, 
rather than hopping parameters t pc i only. As given in 
Table HI upon volume expansion A is reduced by ~ 16%, 
giving rise to a compensating effect of the reduction in 
hopping interaction. This leads to only a reduction by 
less than 4% in covalency energy gain upon volume ex- 
pansion. 

This above expectation is confirmed by our calcula- 
tions on the model Hamiltonian. We carried out two sets 
of calculation, one excluding H e _ e and another including 
H e _ e . In both set of calculations, we assumed Jh = 0, 
therefore Hund's coupling was not operative. The results 
are shown in Fig. [TJ As discussed in the above, the vol- 
ume expansion was found to reduce the covalency energy 
gain due to off-centric movement marginally. Inclusion of 
correlation effects between e g electrons, keeping Jh = 0, 
lead to a further suppression of the volume effect. This 
presumably is rationalized by the fact that inclusion of 
Coulomb interactions would tend to suppress hybridiza- 
tion interaction between oxygens and TM states, thus 
also making the effect of volume variation on hybridiza- 
tion energy gain less effective. 
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FIG. 1: Covalency energy in the 5-band model as a function of 
off-centering shift of Mn cation towards the apical oxygen at 
equilibrium and 2% increased lattice parameters, for H e - e set 
to zero. Energies are shifted in such a way that E cov = at 
u = 0. The negative values indicate a gain in covalency energy 
upon off centric movement. The inset shows the difference 
E cov (2%) — E cov (eq.) with and without H e - e but neglecting 
Hund's coupling. 



In the next step, we investigated the role of Hund's 
coupling between mobile electrons and the localized spins 
on Mn ions. For this purpose, we evaluated the energy 
gain induced by the off-centering distortions at selected 
values of Jh, thus emphasizing its effect on the hybridiza- 
tion mechanism. Wc assumed G-type antiferromagnetic 
ordering for ±2 g spins, which is the magnetic ordering ob- 
served for CaMnC>3. A realistic value for the Hund's cou- 
pling in manganitcs should b o 30 ' 37 between 1 and 2 eV. 
We have varied the Jh value from 0.5 to 2 eV. As pro- 
posed by Khomskii^, we find that the presence of local 
spins on the TM ions and their coupling to the mobile 
electrons through Hund's interaction reduces the energy 
gain due to the off centering movement. As discussed in 
RefH the strong covalent bond between TM e g states 
and neighboring oxygen p states implies the formation 
of a singlet state. This is disfavoured by the strong 
Hund's coupling with the localized t^g electrons on TM 
ion; this interaction can be seen as a local Zeeman field 
due to core spins and effectively acts as a "pair-breaker" 
on the singlet stated. This is demonstrated in the left 
panel of Fig. [2j where we plot the energy difference 
AE cov = [E cov (u; J H ) - E cov (u; J H = 0)] as a function of 
the off-centering distortion u. We find this energy differ- 
ence to be positive, indicating that introduction of Jh re- 
duces the energy gain due to the off-centering movement. 
Increasing Jh reduces more and more the covalency en- 
ergy, making AE cov more positive. Inclusion of electron- 
electron correlation within e g manifold, through inclusion 
of H e - e term, reduces the overall covalency energy due 
to off centering movement, as discussed in the context 
of Fig. [TJ However, the relative energy loss in cova- 
lency energy due to the coupling with the underlying lo- 



cal spins, expressed as E cov (u; Jh = 2)/E cov (u; Jh = 0), 
remains similar between calculations excluding and in- 
cluding H e _ e term, as shown in the inset of left panel of 
Fig. f2J The calculations repeated for the 2% expanded 
lattice constant, exhibit the same trend, as shown in the 
right panel of Fig. f2J This leads us to conclude that the 
effectiveness or role of Hund's coupling to be nearly same 
for the equilibrium and the expanded volume. 

The effectiveness of Hund's coupling, on the other 
hand, is found to depend strongly on the magnetic ar- 
rangements of the underlying ti g spins. The above dis- 
cussed set of calculations have been repeated for the fer- 
romagnetic alignment of the ti g spins. The results are 
summarized in Fig. [3] for the equilibrium volume. Com- 
paring the values of AE cov , we find them to be much 
smaller than those calculated for the G-type antiferro- 
magnetic ordering of t2 g spins and shown in Fig. [2J 
suggesting that Hund's coupling is much less effective. 
AE cov is even found to attain small, negative values for 
Jh between 1 — 1.5eV for moderate values of off center- 
ing movement. Inclusion of -ff e _ e , makes the effect of 
Jh even milder, with much smaller values of AE COV7 as 
shown in the inset of Fig. [3] This finding is in accor- 
dance with the predicted jump in polarization of epitax- 
ially strained SrMnOa, when the system undergoes an 
antifcrromagnetic-ferromagnctic phase transition^ 3 -. 
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FIG. 2: Effect of Hund's coupling on the covalency energy 
at equilibrium (left) and expanded (right) volume as a func- 
tion of the displacement u measured by the energy difference 
AE cov , neglecting other Coulomb interactions. The inset in 
the left panel shows the ratio E cov (u; Jh = 2) / E cov (u; Jh = 
0), measuring the relative energy loss with and without H e - e . 



B. Estimation of elastic stiffness. 

total energy 



Volume effect on 



The analysis, so far, has been restricted to the discus- 
sion of covalency energy gain due to off centering move- 
ment of TM cation. In order to arrive to a complete pic- 
ture of polar instability, one needs to take into account 
the restoring force, which is given by the elastic stiff- 
ness associated to the displacement of TM cation. The 
clastic energy, that emerges from the short-range repul- 
sive interactions between ions, for small displacements, 
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FIG. 3: Effect of Hund's coupling on the covalent energy gain 
at equilibrium volume as a function of the displacement it for 
a ferromagnetic ordering of local ti g spins. The inset shows 
the same including H e - e . 
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FIG. 4: The elastic energy (left panel) and force acting on 
Mn ion (right panel) plotted as a function of off centric dis- 
placement, obtained from DFT calculations. 



is expected to go as ~ ku 2 /2, k being the stiffness con- 
stant. We estimated the stiffness constant from total 
energy calculations carried out on plane wave basis. One 
may identify the covalency energy with the band energy 
contribution to the total energy, and the elastic energy 
as the remainder, i.e. E to t = E cov + E e i as . The relevant 
stiffness constant can then be extracted via a fitting pro- 
cedure from the evolution of E e > as ciS ci function of the 
shift u. This is shown in left panel of Fig. [4] for both the 
equilibrium and expanded volume. The fitting function 
was chosen as f(u) = ku 2 /2 + cu A . As shown in Ta- 
ble UH while the coefficient of the quartic term does not 
change in the expanded volume system with respect to 
the equilibrium case, the stiffness coefficient k is reduced 
by a factor of 0.4 in the expanded volume, suggesting a 
strong reduction of the repulsive short-range interactions. 
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k (eV/k 2 ) 


12.365 


4.961 


10.846 


4.291 


c(eV/A 4 ) 


140 


140 


45 


45 



TABLE II: Stiffness constants as evaluated from estimates 
from elastic energies at equilibrium and 2% increased lattice 
parameters (shown in first and second columns respectively) 
and as evaluated from force acting on Mn cation (shown in 
third and fourth columns respectively). 



In order to check the goodness of our estimate, we also 
adopted an alternative procedure of extracting the stiff- 
ness constant from the force acting on the TM ion. Right 
panel of Fig. 2] shows the evolution of the force acting on 
Mn ion, Fmu as a function of the off centered displace- 
ment. The calculated data points were fitted with the 
function —f'(u) — —ku — 4cu 3 . As shown in Table [TTl 
the estimated k coefficient shows same order of magni- 
tude and same trend under volume expansion, whereas 



the coefficient of the quartic term is more than three 
times smaller when evaluated from the total force. This 
finding, however, is not surprising if we notice that, in 
this second approach, we are neglecting forces acting on 
apical oxygens when we artificially shift the Mn cation. 
We, however, do not expect the physical picture to be 
change much between the two procedures, since in both 
the approaches, while the coefficient k was found to show 
a significant reduction under volume expansion, the co- 
efficient of the quartic term, c, was found to remain un- 
changed by volume expansion. 

From the knowledge of the stiffness constant and co- 
valency energies, in the next step, we investigated the 
ferroelectric instability in cubic perovskite structure of 
CaMnC>3 which arises due to the competition between 
the covalency energy and the elastic energy in the model 
Hamiltonian framework. With the estimated parame- 
ters discussed in previous sections, we find a weak ferro- 
electric instability for the system at equilibrium volume, 
which gets weakened more and more as we allow for a 
finite Hund's coupling between mobile electrons partici- 
pating in the hybridization mechanism and the local Mn 
spins, finally nearly disappearing for Jh = 2 eV , as 
shown in left top panel of Fig. [5j Upon increasing vol- 
ume, however, a strong instability appears, which sur- 
vives even after inclusion of Hund's coupling, though 
gets shifted to lower values of off centric displacement, 
as shown in left bottom panel of Fig. [5] In the above 
calculations, iJ e _ e was set to zero. Inclusion of H e - e , 
as discussed in the above, weakens the covalency effect 
and the ferroelectric instability completely vanishes for 
the equilibrium volume, as shown in the right panels of 
Fig. [5] The ferroelectric instability, however, remains in 
the expanded volume, even for largest Hund's coupling 
strength of 2 eV. 
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FIG. 5: Total energy calculated from the model Hamiltonian 
as a function of the off-centering shift of TM ion at equilib- 
rium (top panels) and 2% increased lattice parameter (bottom 
panels) for different choice of Hund's coupling. Correlations 
within e g electrons have been set to zero in left top and bot- 
tom panels, while they are included at the mean-field level in 
right top and bottom panels. 



IV. CONCLUSIONS 

To conclude, we have introduced a minimal 5-band 
d — p model with Coulomb interactions on the TM ion to 
describe the possible hybridization mechanism for ferro- 
electricity in magnetic Mn-based perovskites. As a case 
study, we have considered the cubic perovskitc structure 
of CaMnOs and have considered two different volumes of 
the structure, the equilibrium volume and 2% expanded 
volume. Our study shows that the ferroelectric instabil- 
ity crucially depends on the balance between two main 
energy contributions, one coming from the tendency of 
empty d states to form covalent bonds with surrounding 
O states and the other from the short-range forces be- 
tween ions, resulting in an elastic energy associated to the 
TM off-centering shift. We found that volume expansion 
has marginal effect on the covalency energy gain, since 
the reduction in hopping interaction gets compensated 
by the simultaneous gain due to change in charge trans- 
fer energy. The coupling of the mobile electrons with 
localized spins on magnetic ions through Hund's interac- 
tion destabilizes the gain in covalency energy due to off 
centric movement, as has been suggested by KhomskiiS.. 
Interestingly, we find the influence of the Hund's coupling 
to crucially depend on the spin arrangements of the lo- 
calized t2 g spins, suggesting a rather strong spin-phonon 
coupling, as predicted in epitaxially strained SrMnO^^. 



Since the hybridization driven mechanism turned out to 
be weakly volume dependent, the effectiveness of Hund's 
coupling in reducing the gain in covalency energy, which 
is a local effect, is found to be weakly affected by volume 
expansion. On the contrary, we found that the stiffness 
of the M n off-centering displacement strongly depends on 
the volume, being more than halved for a 2% increase in 
lattice parameter. The strong reduction of short-range 
forces, combined with an almost unaffected tendency to 
covalency leads to stronger ferroelectric instability in ex- 
panded volume. The correlation between the e g elec- 
trons has the role of suppressing the covalency effect, 
which washes out the weak ferroelectric instability found 
for the equilibrium volume, but retains the much more 
pronounced ferroelectric instability for the expanded vol- 
ume. 

Our results highlight the crucial role played by the bal- 
ance of hopping integrals and charge-transfer energy in 
determining the effectiveness of the covalency. For in- 
stance, much weaker ferroelectric instability predicted for 
the iso-electronic d 3 LaCrOs can be ascribed to a larger 
charge-transfer energy A, as suggested in Ref. HtJ On 
the other hand, one can speculate that the tendency to 
form covalent bonds is reduced for systems with partially 
filled t2 g states, e.g. with d 1 or d 2 , even if Hund's cou- 
pling is expected to play a less effective role compared 
to d configuration; in these systems, the empty t2 g elec- 
trons would establish 7r bonds with neighboring oxygens, 
implying smaller hopping integrals, whereas the charge- 
transfer energy is generally expected to be larger than 
that for d 3 systems. The critical balance between dif- 
ferent energy contributions needs to be estimated as in 
present case. It would, therefore, be interesting to repeat 
the above analysis for other d n perovskite with n < 3. 

Finally, it will be worthwhile to comment on the fol- 
lowing. The ferroelectric instability in perovskite oxides 
is usually explained^ in terms of a close competition be- 
tween short-range (SR) forces, which tend to favor the 
non-polar cubic structure, and the long-range (LR) inter- 
actions, which tend to favor the ferroelectric phase. In 
particular, in Ref. [Usuch SR and LR contributions have 
been discussed in the context of CaMnOs. It may there- 
fore be interesting to connect the various energy contribu- 
tions discussed in the present work to SR and LR forces. 
Since the Hund's coupling as well as the restoring elas- 
tic stiffness tend to destabilize the ferroelectricity, they 
may be linked to SR forces. As detailed in Ref. [2l|, the 
LR contribution has been shown to be directly related to 
anomalously large Born effective charges (see also Ref. 
[39h . We note that the presence of anomalously large 
charges originates from a change in the hybridization be- 
tween occupied and unoccupied electronic states^. In 
the context of CaMnOs, this translates to the change in 
hybridization or covalency between empty Mn c g states 
and occupied O-p states, which is in turn determined by 
both the hopping interaction, t p d, and the charge transfer 
energy, A, both affecting the LR forces. 
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